A cluster mode— coupling approach to weak gelation in attractive colloids 
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Mode-coupling theory (MCT) predicts arrest of colloids in terms of their volume fraction, and the 
range and depth of the interparticle attraction. We discuss how effective values of these parameters 
evolve under cluster aggregation. We argue that weak gelation in colloids can be idealized as a 
two-stage ergodicity breaking: first at short scales (approximated by the bare MCT) and then at 
larger scales (governed by MCT applied to clusters). The competition between arrest and phase 
separation is considered in relation to recent experiments. We predict a long-lived 'semi-ergodic' 
phase of mobile clusters, showing logarithmic relaxation close to the gel line. 

PACS numbers: 



Hard-sphere colloids with short-range attractions can 
undergo several types of arrest. At high densities they 
show two distinct glass transitions (repulsion-driven and 
attraction-driven), with a re-entrant dependence on at- 
traction strength [jj . This scenario was first predicted by 
mode coupling theory (MCT) 0, 0, 0] , and depends on 
both the attraction range <5 (in units of particle diameter) 
and well-depth e (in units of ksT). MCT is remarkably 
successful, at least for large volume fractions tj> > 0.4. 

At lower volume fractions, however, there is no com- 
parable theoretical framework. Yet 'weak gelation', in 
which bonding is strong but not so strong as to be ir- 
reversible, can lead to noncrgodic soft solids, of nonzero 
static elastic modulus, at volume fractions of just a few 
percent It might be argued that a finite modulus 
requires a percolating network of bonds whose lifetime 
exceeds that of the experiment. However this is simplis- 
tic: as shown by the case of repulsive glasses, a finite 
modulus can arise with no bonding at all. We argue here 
that the rigidity of weak gels arises not from bond perco- 
lation but from kinetic ergodicity breaking 6], just as it 
does in glass formation. This suggests that an MCT-like 
approach to weak gelation could be fruitful. 

MCT takes its structural input from equilibrium liq- 
uid state theory; it cannot address states of arrest where 
this structure is strongly perturbed [jj. This matters 
relatively little at <fi > 0.4, where each particle inter- 
acts with many others, and not much room is left for 
structural development upon a quench. But more severe 
consequences must be expected at low volume fractions 
where strongly nonuniform, ramified gels arise. Here the 
pathway to complete nonergodicity (starting from a ho- 
mogenized fluid sample, say) must involve a nontrivial 
episode of structure formation, akin to irreversible cluster 
aggregation (IC A) . Such kinetics certainly dominates for 
irreversible ("strong") gelation (e _1 = 0), where particles 
aggregate on contact into clusters, with various kinetic 
universality classes 0. Relative simplicity is restored at 
low <j) thanks to the invariance of this aggregation process 
under coarse graining in the (ordered) limit 5, e^ 1 — » 
and <j> — > 0. This scaling limit is controlled by an ICA 
'fixed point', where details of the short-ranged attraction 
are irrelevant. The resulting fractal clusters grow indef- 



initely only if (j) = 0; for <fi > they eventually form a 
percolating gel of locally ICA-like structure . 

In this Letter, we explore how ICA might connect 
to weak gelation through a unified (albeit speculative) 
MCT-based scenario of colloidal arrest. We consider the 
changes in the system parameters (0, 5, £ _1 ) upon coarse 
graining in the vicinity of the ICA fixed point. This leads 
to a schematic description of the suspension in terms of 
an effective theory for a dense liquid of 'renormalized par- 
ticles' or coarse-grained clusters, for which recent simula- 
tions [Tol ] seem to provide some direct evidence. Applying 
MCT to this theory, we obtain a new condition for ar- 
rest of the clusters, and identify the arrested state with 
the weak gel phase. If this condition is not met, one has 
instead a fluid of clusters (or 'cluster phase' 11]) that 
is ergodic at large scales. We allow for bond-breaking 
and reconnection at the scale of entire clusters, but not 
reconstruction at shorter scales. We return to this later, 
where we also address the potentially complex interplay 
between weak gelation and phase separation. 

Within our scenario, weak colloidal gelation emerges 
as a double ergodicity breaking: once when the original 
suspension becomes kinetically unstable against aggrega- 
tion, and a second time when the fluid of clusters arrests. 
We suppose that the criterion for the initial instability is 
given by bare MCT @. We assume 1 <C e < oo: bonds 
can break with a small but finite probability. Initially 
this has almost no effect on the aggregation, which flows 
towards the ICA fixed point. But as clusters grow, the 
chances of fragmentation increase. The ICA fixed point 
is unstable and the sytem eventually flows away from it. 
The flow takes place in a parameter space comprising 
effective system variables (t^^e^ 1 ) (Fig.^l. 

These coarse-grained parameters describe the density 
and interaction of aggregated clusters. We write: 



4> = <j) A 3 / dc_1 
5= SN-" 

e = e-x^N- f(N) 



(1) 

(2) 
(3) 



and call the combination of MCT with Eqs.Q]|3J) 'clus- 
ter mode-coupling theory' (CMCT). TV is the number of 
particles in a typical cluster, which is a parametric la- 
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bel along flow trajectories, each of which is indexed by 
values of (<ft, <5, in the initial state after the quench. 
In Eas. 1)112^ 1 . evolution of the effective volume fraction of 
clusters (viewed as quasi-spherical, smooth objects) in- 
volves a scaling index d c , and that of the effective relative 
range S of their attraction another positive exponent v. 

More subtle is the evolution of the bond energy e in 
Eq.©. Since the ICA fixed point leads only to singly- 
connected bonding, there is no power of N multiplying 
the short-range attractive part of the bond energy. How- 
ever, Eq.JJJl accounts for bond breaking through a renor- 
malization of the cluster-cluster attraction e by a break- 
ing entropy, which we write as the logarithm of a breaking 
degeneracy N x < N. This acknowledges that a ramified 
cluster falls apart when any of the N x bonds along its 
backbone is broken (with probability ~ e _£ ). Eq.J3J) also 
includes a term f(N) arising from additional interactions. 
For purely short range attractions this term should van- 
ish, but it allows us to address also the case of a weak 
long range repulsion, such as has recently been shown 
to arise from poorly screened dissociated surface charges 
even for supposedly neutral colloids in organic solvents 
[l2| . Depending on the strength and range of the repul- 
sion, the resulting f(N) may dominate the logarithmic 
term in Eq.@ without having much effect on the equi- 
librium phase diagram p| . 

Initially, all three exponents introduced in the CMCT 
equations are properties of the ICA fixed point (typically 
d c w 2, v m x ~ l/d C ! m three dimensions However, 
they will start evolving once the flow moves away from it, 
as the internal structure of clusters evolves. An explicit 
treatment of this drift is not attempted here: we treat 
d c , v, x as negotiable constants. 

The resulting CMCT flow within the renormalized pa- 
rameter space (i^,(5,£ _1 ) is sketched in Fig. ^ The flow 
is bounded by the MCT transition surface in this space, 
which divides ergodic (upper left) from nonergodic re- 
gions. This surface @>y>3 comprises a dome joined to 
a sheet. Below the dome lie attractive glasses and to the 
right of the sheet, repulsive ones. Note the presence of a 
seam <j) c , on the transition surface, along the line where 
the sheet meets the dome. Consider a trajectory of in- 
creasing N which rises up through the dome, left of the 
scam (cj) < 4> c ). At this point, aggregation ceases, be- 
cause the clusters' effective parameter values now corre- 
spond to an ergodic phase. A semi-ergodic fluid of finite 
clusters can be expected whenever the initial parameter 
values belong to such a trajectory. (This requires small 
8.) If a series of such cluster fluids is now created by 
increasing <f> at fixed e (or vice versa), then the last of 
these to enter an ergodic phase will do so just beside the 
seam (<p — > <f>~). Accordingly the first fully nonergodic 
phase formed has a trajectory that just meets this seam; 
beyond lies a repulsive glass, so that aggregation ceases. 
Unlike the cluster fluid, the resulting marginal gel has 
finite modulus — as do the non-marginal gels arising at 
larger (f> or e. This not because of percolation, which 
for reversible bonding is not relevant, but because the 
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FIG. 1: Schematic parameter flows (arrows) near the ICA 
fixed point (at the origin, bottom left). Flow (a) within the 
lower coordinate plane (e _1 = 0) gives percolation of ICA 
clusters on passing through a non-equilibrium percolation line 
(heavy dashed). For f(N) > 0, trajectories starting from 
finite e~ flow upwards to smaller e and towards larger <f> and 
smaller 5. Those that arrive on the domed dark surface result 
in semi-ergodic cluster phases. The marginal gel trajectory 
(b) ends at the seam where the repulsive (light surface) and 
attractive (dark surface) branches of the MCT transtion meet, 
close to a line of A3 singularities (dash-dotted) or/and an 
A4 singularity (*). Inset: <j> — i~ —projection of trajectories 
leading to cluster or gel phase (thin) and marginal gel (b). 



CMCT flow has led to a "cluster glass" which we iden- 
tify with the weak gel phase. A typical phase diagram 
predicted by this approach is shown in Fig. [5J 

Note that for moderate 5, the seam is not too far from 
the line of ^3 singularities, which is associated with a dis- 
tinctive arrest scenario including a large exponent for the 
divergence of the structural relaxation time, and strongly 
stretched exponential or even logarithmic decay of the 
density correlations 0, Q, 0] . Marginal gels could thus 
be influenced by these singularities even at low volume 
fractions, although the details of this depend on the ini- 
tial range 5 and on the exponent v. 

CMCT clearly offers nontrivial predictions for cluster 
phases, weak gels and relaxation anomalies. But it ne- 
glects two important pieces of physics which may com- 
plicate the experimental situation considerably. Firstly, 
we allowed clusters to dissociate by bond breaking, but 
neglected any internal restructuring within the clusters. 
By allowing these to coarsen into globules, this could 
cause both the cluster fluid and weak gel phases to have 
finite lifetime. Moreover, for f(N) — the breaking 
and rebonding of the O(N) dangling ends within a clus- 
ter may cause local restructuring even before the flow 
reaches the MCT surface. This issue is less important in 
the presence of a weak, long-range repulsion f(N) > 
(e.g. |12|). which, like bond-breaking, drives an upward 
flow in Fig. ^ In favorable cases this could dominate 
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FIG. 2: Kinetic phase diagram predicted from CMCT with 
initial range of square-well attraction S = 0.1. Exponents are 
chosen for definiteness as d c — 2, \ = v — 1/2 and f(N) = 
0. Bare MCT (in the low-<f> approximation 6]) is used to 
estimate the onset of a clustering instability (lower solid line) . 
The CMCT line divides cluster fluids from weak gels (upper 
solid line). Results allowing also for phase separation, using 
the indicative model of Fig. |3] give shifted phase boundaries 
(dashed lines). (The part of the equilibrium binodal above 
the AB tieline in Fig. |3]is also shown, dotted.) 



the bond-breaking term at large e, and allow the CMCT 
scenario to establish on a timescale during which local 
reconstruction remains negligible. 

The second, closely related, omission from the CMCT 
predictions of Fig. [3 is the interplay of cluster forma- 
tion with liquid-gas type phase separation 0. To be 
confident of observing the CMCT line directly at small 
<j>, it must lie outside the binodal. When that condi- 
tion is not satisfied, we expect results that depend on 
quench rates. For a very slow quench, the system first 
demixes locally and then coarsens, until the denser coex- 
isting phase, whose properties continue to evolve, arrests 
by a CMCT-like mechanism (merging with bare MCT at 
high enough density). Slow quenches have recently been 
examined experimentally, and these ideas may account 
for a surprising "bead phase" seen at low <f> [F7|. 

For faster quench rates, local arrest can precede coars- 
ening so that at small <j>, large fractal clusters will at first 
evolve, regardless of the underlying phase behavior. The 
resulting scenario is explored, within an indicative model, 
in Fig. |3| We follow Keller [ljj in arguing that the dens- 
est phase capable of formation by demixing arises where 
the glass line intersects the binodal (point B). The non- 
ergodicity line should depart from the CMCT prediction 
at this point, and track to lower densities along the tie 
line BA. This describes a locally demixed state in which 
only the denser (gel-like) B-componcnt is nonergodic. 
At high <fi the B-phase percolates giving a finite modu- 
lus gel; but at lower 4> one expects a cluster fluid, now 
made of disconnected gel-like fragments. The 'gel line', 
separating states of finite modulus from those without, 




FIG. 3: Kinetic phase diagram, superimposed on the equi- 
librium phase diagram, for an indicative model. The equi- 
librium free energy density is that of a lattice gas, f(4>) = 
<f>\n(f> + {4> m — 4>) ln(<^7n — 4>) — \<j>{<j>m — 4>) with <f) m = 0.64 
chosen as a maximum packing fraction. The second virial is 
matched to the colloidal system so that A(e) = 1/2 — B2(e)v 
with v the volume of a colloid; B2 is that of a square well 
potential of range 5 = 0.1. The CMCT line is found for the 
colloid system as in Fig. [5] equilibrium binodal (Bi) and spin- 
odal (Sp) are marked. Inset: modified tangent construction 
on f(cj>) |Ts| in which states beyond D are deemed inacces- 
sible. The gel line (with percolation threshold set at phase 
volume 0.3 of the arrested phase) bounds the darker shaded 
region. Between this and the modified binodal (dashed at 
left) lies the cluster fluid phase (lighter shading) . 

corresponds to the geometric percolation threshold of the 
arrested phase. Deeper quenches should behave similarly, 
with local coexistence between a dilute fluid C and an 
arrested state D whose density lies on the CMCT line. 
This is governed by a modified common tangent |18| in 
which inaccessible densities (beyond point D) are prohib- 
ited (Fig. El inset) shifting the binodals inward. Along 
the tieline CD, the gel line again corresponds to percola- 
tion of the arrested phase. At still deeper quenches the 
CMCT line lies within a (renormalized) spinodal curve; 
this local instability complicates matters further, and is 
ignored in the indicative model of Fig. [21 But a transi- 
tion from percolating gel to a cluster phase of gel frag- 
ments should remain. Finally, at infinite e the CMCT 
line, the gel line, the equilibrium spinodal, and the bin- 
odal all meet the ICA fixed point at the origin. Notably, 
this scenario allows the 'marginal gel' behavior to extend 
inside the gel phase: indeed a range of compositions, in- 
cluding the AB tieline and a region below it, now have 
local phase separation involving a marginal gel (Fig. ^ . 

We now confront the implications of our work with 
experimental data. The most direct consequence of our 
hypothesis of a double ergodicity breaking, a fractal gel 
structure, is a very common observation. It has inspired 
successful phenomenological gel models, e.g. [l9l I20L l21| . 
for which CMCT could provide a more rational (albeit 
still speculative) basis. Experiments demonstrate 
that even for moderate attractions the assumed closeness 
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to the ICA limit is well justified at low volume fractions 
in some cases, while in others the measured fractal di- 
mension lies higher than for the appropriate ICA model 
[23|. The latter could well be a consequence of a partial 
phase separation followed by arrest, as discussed above. 

Secondly, CMCT predicts dynamic anomalies along 
and near much of the gel line (including the part that lies 
along AB tieline of Fig. [3J. In this context, the observa- 
tion in Ref.^l of a large exponent (5.5±1 at (f> « 0.1) for 
the divergence of the terminal time, and of slow, possibly 
logarithmic, structural relaxation close to the gel line at 
low (f> are both suggestive of CMCT. 

A final comparison involves the nonequilibrium phase 
diagr am, where a striking feature at low concentrations 
(e.g. ^3 0, 0|) is indeed the existence of a phase of 
large mobile aggregates, lying between the ordinary fluid 
phase and the gel. In fact, the phase diagrams in Fig. 5 
of Ref . an d Fig. 2 of Ref . ||2f| are very similar to our 
Figs. |3 El respectively. As mentioned above, an inter- 
play of CMCT with phase separation is expected for 
most colloidal parameters, including those of Ref. [ill] . 
Fig. [21 therefore shows also the phase boundaries ad- 
justed to allow for phase separation according to the 
indicative model of Fig. |3| In this representation the 
cluster/gel boundary remains nearly straight, although 
shifted downward somewhat. The results remain qualita- 
tively consistent with Ref.^J; given the indicative char- 
acter of the model, this is encouraging. Since the system 



is among those for which a long range Coulomb repulsion 
might be present 0], it is possible that f(N) > 0. Qual- 
itatively, this would move the calculated gel line upward, 
stabilizing the cluster phase. 

In summary, we have proposed a new framework to 
describe weak gelation based on mode-coupling theory 
applied to clusters (CMCT). This certainly contains spec- 
ulative elements but gives nontrivial predictions, several 
of which compare well to existing experimental data. Un- 
derlying our picture of a double ergodicity breaking is the 
premise of a time scale separation. Our "cluster glass" 
view of weak gels is a meaningful concept on time scales 
longer than the time for cluster formation, but short com- 
pared to typical times over which a cluster reconstructs 
as it evolves towards equilibrium. This condition is more 
easily satisfied in systems where a repulsion between clus- 
ters builds up as they grow (e.g. weakly charged parti- 
cles), limiting their final size. When the size is limited by 
bond-breaking, the required time scale separation is less 
likely. But even if the resulting weak gels are not true 
examples of broken ergodicity and are only temporary 
structures, CMCT may offer valuable insights into their 
behavior. Further experimentally relevant consequences 
of the scheme will be the subject of future work. 
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